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GCIRS16SW: A MASSIVE ECLIPSING BINARY IN THE GALACTIC CENTER * 

F. Martins \ S. Trippe \ T. Paumard \ T. Ott \ R. Genzel L2 , G. Rauw 3 , F. Eisenhauer 1 , S. Gillessen \ H. 

Maness 4 , R. Abuter 5 

ABSTRACT 

We report on the spectroscopic monitoring of GCIRS16SW, an Ofpe/WN9 star and LBV candidate 
in the central parsec of the Galaxy. SINFONI observations show strong daily spectroscopic changes 
in the K band. Radial velocities are derived from the He I 2.112/im line complex and vary regularly 
with a period of 19.45 days, indicating that the star is most likely an eclipsing binary. Under various 
assumptions, we are able to derive a mass of ~ 50 M© for each component. 
Subject headings: Stars: binaries: eclipsing — Stars: early-type — Galaxy: center 



1. INTRODUCTION 

The central cluster constitutes one of th e largest con- 
centr ations of massive stars in the Galaxy ijGenzel et alJ 
|2003|) . Nearly 100 OB and Wolf-Rayet stars are confined 
in a compact region of radius ~ 0.5 parsec centered on the 
super- massive black hole asso ciated with the radio source 
SgrA* jPaumard et a l. 2006). Among this population of 
young massive stars, six are thought to be Luminous 
Blue Variables (LBV): IRS16NE IRS16C, IRS16NW 
IRS16SW, IR S33E and IRS34W llPaumard et alJl200l 
iTrippe etaDl2006|) . LBVs are evolved massive stars ex- 
periencing strong variability in both photometry and 
spectroscopy du e to their proximity to the Hu mphreys- 
Davidson limit {Humphreys fc Davidsonl[l9 94) , a region 
of the HR diagram where the luminosity of the stars 
reaches the Eddington luminosity so that instabilities 
develop in their atmospheres, leading to strong mass 
ejection and drastic changes in the stellar properties 
(T e s, radius). The six stars mentioned above are only 
LBV "candidates" (LBVc) since they have not been ob- 
served to experience the strong outbursts and photomet- 
ric changes typical of bona fi de LBVs such as 77 Car 
{Davidson &: Humphreys! 11997^1 . However, their lumi- 
nosities and spectra are very similar to stars known to 
be "quiescent" LBV, i.e. stars having experienced an 
LBV event in the past and being now in a more sta- 
ble phase. In addition, one of them - IRS34W - has 
shown photometric variability on timescales of months 
to years which was interpreted as the formation of dust 
from material prey iously ejected by an LBV outburst 
{Trippe et al]l2006|) . 

Among these six stars, IRS16SW deserves special at- 
tention. This star was c l aimed to be a massive eclips- 
ing binary by lOtt et akl {1999ft since its K band mag- 
nitude displays regular variations with a periodicity of 
9.72 days. However, the absence of a second eclipse in 
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the light-curve lead iDePov et alJ {2004|) to the conclu- 
sion that the binary scenario was not correct, and that 
IRS16SW was instead a pulsating massive star, a class 
of star predicted by theory but not observed so far. 

Here, we present results of the spectroscopic monitor- 
ing of IRS16SW revealing periodic variations of radial 
velocities which are interpreted as the signature of a mas- 
sive spectroscopic and eclipsing binary . 

2. OBSERVATIONS AND DATA REDUCTION 

We used SINFONI (lEisenhauer et alJ 12003]) on the 
ESO/VLT to obtain spectra of IRS16SW. Observations 
were carried out under seeing limited conditions and were 
performed on August 28 th - September 1 st , September 
4 th , October 2 nd , and October 4 th - 12 th 2005, and March 
18 th - 21 st 2006. In order to get the best spectral reso- 
lution available with SINFONI (R = 4000) we restricted 
ourselves to the K band. Short exposures (2 x 60 seconds) 
were sufficient to obtain S/N ~ 30. Da ta reduction was 
performed as in lEisenhauer et all (|2005f) . The final spec- 
tra were subsequently carefully extracted from the "data 
cubes" by selecting a circular aperture (radius of 3 pixels) 
centered on the star and by subtracting from it an annu- 
lus of inner (outer) radius 3 (4) pixels. This procedure 
allowed a good removal of nebular contamination. 

3. RESULTS 

3.1. Spectroscopic variability and radial velocities 

Fig. n shows the variation of the He I 2.112 \im and 
Br7 lines with time. It is obvious that not only the line 
shape but also the position of the centroid varies. In or- 
der to test the binary scenario, we have derived radial 
velocities (RV). For that purpose, we have used the line 
at 2.112 fim since it is less affected by wind emission 
than other lines and is formed closer to the photosphere, 
allowing a better estimate of the star's motion. This line 
is however a blend of at least two He lines, and synthetic 
spectra computed with atmosphere models reveal that 
the position of the strongest absorption part of the pro- 
file can vary by several 100 km s _1 around 2.112 \xra 
depending on the stellar and wind properties. We have 
thus adopted this wavelength as our reference, but we 
stress that the absolute value of the RV may be system- 
atically shifted compared to the real value due to this 
choice. In practice, we have measured the position of 
the maximum absorption trough in the line complex and 
computed the radial velocity from the wavelength shift 
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Fig. 1. — Montage of spectra of IRS16SW around the region of He I 2.112 fim (left) and Br7 (right) taken between October 2 nd and 
October 12 th 2005. The chan ges in both the line shape and position are clearly seen. 



compared to the adopted reference wavelength. When 
the line shows a double peak, we have always measured 
the position of the deepest absorption part of the profile 
(which also turned out to always be the bluest). This 
implicitely assumes that if the star is a binary, this ab- 
sorption is always produced by the same star. Note that 
when present, the two absorption peaks are separated 
by the theoretical spacing between the two Helium lines 
around 2.112 /j,m: 170 km s — 1 . Hence we conclude that 
the second peak most likely comes from the same star. 
The derived RV are presented in Fig. [21 and follow very 
nicely an almost exact cosine curve (see below). This is 
a strong indication that IRS16SW is a binary star (and 
justifies a posteriori our method to derive RVs). 

A period folding analysis applied to the RV curve gives 
19.3 ± 0.4 days, in good agreement wit h the K band 
light-curve analysis. iDePov et al.1 l|2004f ) derived a pe- 
riod of 9.725 ±0.005 days but argue that if this light- 
curve was to be produced by a binary star, the absence 
of second minimum should point to a system composed 
of two stars with the same K band luminosity, and con- 
sequently to a true period of 2 x 9.7 25 = 19.45 days. A 
re-analysis of the photometric data of lOtt et alJ 1^)99) in 
view of these new results confirms that a period of 19.447 
±0.011 days is indeed present in the period folding di- 
agram (for a description of the method, see lOtt et all 
1999). We conclude that IRS16SW is most likely a sin- 
gle line spectroscopic (SB1) and eclipsing binary with an 
orbital period of 19.45 days. 



3.2. Orbital solution and physical parameters 

With the RVs in hand, we have performed an or- 
bital solution for an SB1 binary using the method of 
iRauw et al.l l)2000fl . which is based on t he Wolfe, Horak 
& Storer algorithm ijWolfe et al.lll967|) . The resulting 
parameters are given in Tabled Note the small eccen- 
tricity justifying that the fit of the RV curve in Fig. is 
almost a cosine. The orbital solution also gives the mass 
function 



f(m) = 



(M 2 sini) 3 
(Mi + M 2 ) 



10.47M, 



(1) 



from which one can estimate the individual masses of 
each component (Mi and M 2 ) if one uses in addition 
Kepler's third law 



2^3 



Mi 



Mo = 



in 2 r 
GP 2 



(2) 



where r is the separation between the two stars and P 
the period. However, one needs an estimate of 1) the in- 
clination and 2) the separation. As for the latter, the ab- 
sence of plateau in the K band light-curve indicates that 
contact is achieved in the binary: as soon as the primary 
eclipse ends, the secondary eclipse starts. Hence, one can 
assume that r is simply the sum of the radii of the two 
components. The value of asini we find (and a itself, 
sun being close to 1, see below) i s similar to the stella r 
radius of LBVc stars derived by iNaiarro et al.1 fl997), 
so we assume r = 2 x a. In that configuration, the two 
stars, which we know have similar K magnitude, have 
similar radii and rotate around each other, the center 
of mass being the contact point. To get an estimate of 
the inclination, we have used the software NIGHTFALL 
6 to fit the light-curve (see result in Fig. EJ. Assum- 
ing a similar effe ctive temperature of 28000 K for both 
components (see INaiarro et al.l 119971) and a light ratio 
of one, we obtained a reasonable fit with an inclination 
i ~ 70 deg. Note that we had to adopt Roche lobe filling 
factors of 1.3 (the maximum allowed value in NIGHT- 
FALL) to correctly reproduce the light-curve. If filling 
factors lower than one are used, the light-curve can not 
be reproduced correctly as the "peaks" are too broad (a 
problem encountered bv lDePov et al]l2004D . This con- 
firms that contact is achieved and justifies our assump- 
tion that the separation is the sum of the stellar radii. 

' software developed by R. Wichmann 
and freely available at the following URL 
http: / /www. hs. uni-hamburg.de/DE/Ins/Per/Wichmann/Nightfall. html 
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TABLE 1 

Orbital parameters - semi-major axis, eccentricity, 
systemic velocity, amplitude, longitude of periastron - as 
derived from the analysis of the radial velocity curve. 



Fig. 2. — Radial velocity curve of IRS16SW together with the 
best orbital solution (for P = 19.45 days, solid line). The typical 
uncertainty on the radial velocity is ± 20 km s _1 . Parameters for 
the best fit solution are given in Table 111 




1 

phase 



Fig. 3. — K band light-curve displayed for a period of 1 9.45 days, 
togeth er with the best fit (dashed line). Data are from lQ*tt et al.l 
11999). The two eclipses are similar, indicating that both compo- 
nents of the binary have the same K band magnitude. Note also 
the absence of plateau which suggests that contact is achieved. 



Given the limitations on the light-curve fit and the 
complications in the physics of the star due to contact 
(mass transfer, departure from sphericity due to gravi- 
tational interaction, hot spot in interaction region), we 
stress that our estimate of the inclination is only indica- 
tive. It is however consistent with the value expected for 
an eclipsing binary. With i ~ 70 deg, and using Eq. ^ 
and |3 we finally derive Mi ~ M 2 ~ 50M Q . 

3.3. Spectral disentangling 

In order to get more insight into the properties of the 
components of IRS16SW, we have att empted to disentan- 
gle th e spectra using the method of IGonzalez fc Levatol 
(2006). In practice, the RVs of the primary are used to 
evaluate an average spectrum in the primary's rest frame. 
This provides an approximation of the primary spectrum 



asini [Rq] 


66.4 ± 2.2 


e 


0.088 ± 0.023 


vo [km s _1 ] 


459.5 ± 3.6 


Kl [km s" 1 ] 


173.8 ± 5.5 


w [deg] 


334.0 ± 18.3 




Fig. 4. — Spectra of the primary (top) and secondary (bottom) 
components of IRS16SW as obtained from the spectral disentan- 
gling analysis. These spectra should be interpreted with caution 
since they were obtained under several assumptions (see text). 



which is then shifted back into the observer's frame and 
subtracted from the observed spectrum. The residual is 
composed of the secondary spectrum from which RVs can 
be estimated. The whole procedure is then re-started in- 
verting the role of the two components, and is iterated 
until convergence. 

Due to the limited number of spectral lines and S/N 
ratio of our spectra, no solution could be found leaving 
all parameters free. We therefore decided to freeze the 
primary RVs at their values determined above. Given 
the results of the orbital solution, we also decided to set 
the mass ratio to 1.0, and assumed that both compo- 
nents have the same systemic velocities. In that case, 
convergence could be achieved and the resulting spectra 
are shown in Fig. These spectra should be interpreted 
with caution. Not only they were obtained under the 
assumption that the mass ratio is 1.0, but the procedure 
used also implies that the spectra are free of contamina- 
tion by wind-wind collision or any other interaction in 
the contact region. With these restrictions in mind, the 
main qualitative conclusion we draw is that the spectra 
of both stars, and consequently their properties, are sim- 
ilar and typical of LBV candidates such as the Pistol star 
l|Figer et alJfl999[) . 

4. DISCUSSION 
4.1. Binary versus pulsating variable 
iDePov et all l)2004|) argue that IRS16SW was a pulsat- 
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ing massive star based on the absence of second minimum 
in the light-curve and the difficulty to fit this light-curve 
in the binary scenario (but see Sect. 13.21) . They com- 
pare the observed variation in K mag nitude to the pre- 
diction s of the dynamical models of iDorfi fc Gautschvl 
(2000) and conclude that there is a reasonable qualita- 
tive agreement. However, there are some quantitative 
discrepancies. First, the amplitude of the variation is 
much larger than predicted: although IDorfi fc Gautschvl 
( 2000) do not compute K band photometry, one can es- 
timate the variation in this band to be at most 0.2 mag 
(inspection of their Table 2 reveals that the amplitude 
of photometric variations decreases with wavelength and 
is < 0.2 mag in the 1 band), while we observe 0.55 mag. 
Second, the period we derive — 19.45 days — is larger 
than expected in the p ulsating scenario (see Table 1 of 
IDorfi fc Gautsch\ll2(i0Ch . 

Concerning spectroscopic changes, although in princi- 
ple one can not completely rule out the possibility that 
they are due to motions of the atmosphere and fluctu- 
ations of t he physical parameters (Teffi radius) due to 
pulsations ijDorfi fc Gautschvl 12000). the timescales are 
again not consistent: ~ 1 day for pulsations compared to 
19.45 days observed. Besides, so far there are no theoret- 
ical predictions of spectroscopic changes caused by pul- 
sations to which we could compare our observed spectra. 
The binary nature of IRS16SW is thus strongly favored. 

The absence of secondary eclipse in the light-curve is 
explained by the similar K band magnitude of the two 
components. T his is another indica tion that both stars 
are very similar. iDePov et alJ 1)20041) report the presence 
of a variation in H — K on a period of 9.725 days, H — K 
being bluer when the system is brighter. A similar trend 
was observed in the optical ph otometry of the m assive 
contact binary V606 Centauri ijLorenz eFaI]|l99<|. We 
interpret this as a sign of heating in the contact zone, 
making the spectral energy distribution in this region 
bluer. Massive binaries are indeed known to produce X- 
rays through colliding winds, and we might expect the 
same kind of interaction and heating around the contact 
region. Again, since this region is seen twice during an 
orbital revolution, an observed period half the true or- 
bital period is naturally derived from the H — K curve. 

4.2. Stellar evolution and the LBV phenomenon 

Whether or not all massive stars go throu g h the LBV 
phase is still u nder debate. lLanger et alJ l)1994|) and 
iPasauali et alJ lll997|) argue that this is t he case, while 
othe r observational iCrowther et al.lll995|) and theoret- 
ical ijMevnet fc Maederll2005|) studies indicate that the 
most massive stars (M > 60M Q ) may skip this phase. 



This is an important issue since although short, the LBV 
phase is crucial in the mass loss history, and consequently 
in the su bsequent evolutiom of m assive stars. Recent 
studies bv lSmith fc Owockil l)2006|) even claim that most 
of the mass of hot stars is lost during the LBV phase. 
Here, we provide an accurate measurement of the present 
mass of a candidate LBV, confirming that a star with an 
initial mass larger than 50 M© (and likely of the order 60- 
70 Mq) may become a LBV. This is a strong constraint 
for evolutionary models. 

Of course, one could argue that the star's evolution was 
affected by binarity. However, inspection of the spectral 
morp hology and physical properties (see iNaiarro et alJ 
{m% of IRS16SW and the other LBVc in the Galactic 
Center shows similarities. Our monitoring of IRS16SW 
also includes the other LBVc in the Galactic Center. Ex- 
cept for IRS16NE, none of these stars showed any spec- 
troscopic variation, ruling out the possibility that they 
are close bina ries. IRS16NE sho wed some RV fluctua- 
tions (see also lTanner et al1l2006j) . but so far we do not 
have enough data point to sample an hypothetical RV 
curve. Hence, the similarity between the spectrum of 
IRS16SW and those of the other single LBV candidates 
leads us to the conclusion that binarity has not (yet?) 
significantly affected the evolution of IRS16SW. Since 
the mass of the two components is similar, one can spec- 
ulate that IRS16SW is composed of two stars initially 
equally massive that have so far evolved in parallel in a 
detached system, without influencing each other. They 
may have just entered the LBV phase during which con- 
tact was achieved due to their respective expansion. This 
event probably happened very recently (the LBV phase 
lasting ~ 10 5 yr) so that the general properties of both 
components have not yet been affected by mass transfer 
and binary evolution. Such a scenario is consistent with 
both stars displaying similar spectra (but see Sect. 13.31 
for caution words). 

Assuming that 1RS16SW was not affected by bi- 
nary evolution, the properties of the LBV candidate in 
IRS16SW can thus be used to constrain evolutionary 
models of single massive stars. In their rec ent models in- 
cluding rotation. iMevnet fc Maederl l)2005|) stressed that 
the LBV phase is not systematically reached above 45 
M Q . Here, we have an example for which it is the case 
(under the assumption that IRS16SW will turn — or has 
already turned in the past — into a genuine LBV). 
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